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We present a description of nuclear spontaneous fission, and generally of quantum tunneling, in 
terms of instantons - periodic imaginary-time solutions to time-dependent mean-field equations - 
that allows for a comparison with more familiar and used generator coordinate (GCM) and adiabatic 
time- dependent Hartree-Fock (ATDHF) methods. It is shown that the action functional whose value 
for the instanton is the quasiclassical estimate of the decay exponent fulfils the minimum principle 
when additional constraints are imposed on trial fission paths. In analogy with mechanics, these 
are conditions of energy conservation and the velocity-momentum relations. In the adiabatic limit 
the instanton method reduces to the time-odd ATDHF equation, with collective mass including 
the time-odd Thouless-Valatin term, while the GCM mass completely ignores velocity-momentum 
relations. This implies that GCM inertia generally overestimates instanton-related decay rate. The 
very existence of the minimum principle offers a hope for a variational search for instantons. After 
the inclusion of pairing, the instanton equations and the variational principle can be expressed in 
terms of the imaginary-time-dependent Hartree-Fock-Bogolyubov (TDHFB) theory. The adiabatic 
limit of this theory reproduces ATDHFB inertia. 



I. INTRODUCTION 



Decay of a metastable state of a system of interacting fermions or bosons is an important phenomenon relevant 
to nuclear, atomic and condensed matter physics. The calculation of decay rate requires the exact knowledge of the 
wave function in the proper asymptotic region which is usually very difficult to achieve for many-body systems. In 
fact, very often the only feasible description of systems including hundreds or more particles relies on the quantum 
mean-field theory. Unfortunately, such theory does not contain quantum tunneling. This gives rise to a notorious 
arbitrariness in calculations of decay rates or half-lives which concerns a selection of relevant degrees of freedom and 
prescriptions for potential and inertia parameters. 

Specifically, within the Hartree-Fock (HF) method, static equations give only saddle points on energy landscape 

n[r,i']^ f dxj2^'^rk^i^k + nr,i^], (i) 

J A; ™ 

with V[ip* , ip] being potential energy, so one has to resort to the time-dependent HF (TDHF) equations for dynamics 

thdtMt) = h{t)Mt) - -^^'Mt) + (2) 

with the mean-field single-particle (s.p.) Hamiltonian h{t) given by h^ip* {t),'(p{t)]ipk{t) = ^T^/H-'ki^)^ f^nd the self- 
consistent s.p. potential V{t) given by SV/SipKt) — V{t)ipk{t). For the case of energy H given by a density functional 
we assume in the following (if not indicated otherwise) that it has properties of the expectation value of the Hamil- 
tonian. Although Eqs.([2]) look like the Schrodinger equations, in fact, they are classical field equations, due to a 
nonlinear dependence of h on ipk- Energy ([T]) and overlaps {ipk \ ^re conserved by Eqs.Q. The former forbids 
a tunneling within TDHF, i.e. an escape from a minimum of 7i with energy lower than the saddle. Evidently, this 
comes about by projection of the full many-body theory onto Slater states. 

A quasiclassical treatment of quantum tunneling within the many-body mean-field theory, which is a natural 
generalization of the Gamow treatment of alpha-decay to an infinite-dimensional system of fields, leads to instantons 
- periodic imaginary-time solutions to TDHF equations [l|. This method exploits an idea of trajectories evolving 
in imaginary time [4] which emerge from the stationary-phase approximation to the path-integral expression for 
Tr{E — H)^^ . The decay rate of a metastable state is proportional to exp{—S/fi), where S is action for the optimal 
instanton. We do not consider here a prefactor coming from quantum fluctuations around the optimal path. 

For a particle in an external potential such optimal decay trajectory describes classical motion in the inverted 
potential. It starts at the metastable state (being a local maximum of the inverted well) and returns there after 
bouncing from the inverted barrier; hence the name 'bounce'. For a system of interacting fermions, one has to 
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transform TDHF Eqs.(l2]) to imaginary time, i.e^, formally, t —ir. Under this transformation, ■0 tfj{x,—iT) = 
4i{x,t) and ip* —f ^p{x,—iT)* ~ (j){x,—T)* [H, y|- It follows that density p{x,t) = tp* {x,t)^p{x,t) transforms to 
pix,T) = (j){x,—T)*(j){x,T). This has important consequences. First, the mean-field equations in imaginary time 
lii: 

n'-^ir) . - (Hr) - Mr) = |^VV.(r) - ^ + e.Mr), (3) 

are non-local in r, as V as well as /i(t) = (— t), 0(r)] depend on both 0(a;,r) and (j){x,—T). Second, density 
p{x,t), generally complex or piecewise negative, does not correspond to any Slater determinant, unlike in the real- 
time dynamics. In analogy with TDHF, Eqs.([3l) conserve energy H(r) = r), </)(t)], with: 

n{T) = j dxJ2^^<PU-r)VMr)+V[rhr),cP{T)]. (4) 

k 

The above formula means that one obtains Ti-lr) replacing everywhere ipl{t) by t) in the usual form of the 
energy functional. Since the Hamiltonian is hermitean, = H, it follows that 7i(— r) = 7Y*(r) and the mean-field 
Hamiltonian /i(r), defined by h{T)(t>{T) = 6Ti./S(p*{—T), fulfils the condition h{—T) = h^{T). The latter ensures that 
Eqs. ^ without the ek(j)k term conserve the overlaps: -^{4>i{—T) \ 4>j{T)) = 0. The complete Eqs.® still conserve 
diagonal overlaps, while giving the exponential time-dependence to the off-diagonal ones. However, those overlaps 
remain zero for all r, if equal zero at some r. 

As usual, the saddle point approximation to the path integral leads to the periodicity condition for the optimal 
trajectories. Hence, bounce is a periodic intanton, 

MT/2) = M-T/2), (5) 

and the periodicity is enforced by the ek<t>k term in the Eqs. ([3]). The physical context imposes the specific boundary 
conditions on bounce. For a description of the decay of a metastable ground state, the initial (and thus also the final) 
states have to be chosen equal to the HF solutions "0^^ at the metastable minimum, (f>k{T/2) — T/2) = i^k^ , 
with total energy Egg, and the parameters must be equal to the HF s.p. energies at this minimum. The s.p. 
states (f>k{T — 0) form some normal (as (pK—r) = (pKr) at t = 0) HF state at energy H ~ Egg on "the other side 
of the barrier" . The periodicity condition together with the initial condition fix the particular constant values of the 
overlaps: 

(4(-r)|0,(r))=J,,. (6) 

Decay exponent is given by [3, [1| : 

rT/2 

-^/2 k 



S^hf^'^drY^iM-r) I ^(r)). (7) 



Bounce penetrates the static barrier, impermeable for real-time solutions at the same energy, practically in a finite 
time interval around r = and becomes infinitely slow close to the endpoints, so that T extends to infinity P, S 0]- 
Eqs.© dete rmine both decay channels and decay probabilities. No additional assumptions are necessary, as they 
form a complete quasiclassical solution to the tunneling problem within the mean-field theory. 

Up to now, solutions of Eqs.® have been obtained only for relatively simple systems [11,11, 0, S S [3l- The task 
of finding instantons seems rather hopeless without a special treatment: to handle non- locality in r one could try to 
solve Eqs. ([3]) together with 

.n^M_lA + (^ki-r)-e,)M-r)=0, (8) 

describing instantons evolving backwards, obtained from (jSj by using the identity {dTf){—T) = —drifi—T)). However, 
Eqs.® describe the inverse diffusion [cf signs of time and spatial derivatives in ([3]) and ([5])], which leaves no hope 
for a stable solution. The problem seems more difficult than a search for periodic solutions of the real-time TDHF 
equations, which is known to be difficult enough. In the presented form, instanton approach did not lead to any 
comparisons with the existent studies of fission, which are mostly based either on the generator coordinate method 
(GCM) in the gaussian overlap approximation (GOA) or on the adiabatic TDHF (ATDHF) method, either in its 
extreme cranking or some more refined version. 
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In this work we present the instanton method in famihar terms of the HF or HFB theory. This helps to grasp 
similarities and differences between this and other methods and to clarify their interrelations. In particular, the 
instanton turns out to be an analogue of the self-consistent TDHF in the representation of the time-even and time- 
odd components of the density matrix Q. Both produce the same inertia in the ATDHF limit, when one includes 
time-odd components only to the first order. 

Moreover, it turns out that instanton action Eq.l?]) is the minimum value of the action functional over properly 
constrained set of trial fission trajectories defined in the space of Slater determinants Thus, any fission path 
which satisfies these constraints provides the upper bound for the decay exponent. This offers a hope for a variational 
approach to finding instantons. One may also expect that a good estimate of action may be easier to find than that 
of the instanton itself. 

The starting point is the realization that Eqs.® describe two different sets of Slater determinants, bra t) 
built out of r), and ket ^{t) built out of (j)k{T)j while energy Ti. is equal to the off-diagonal energy overlap 
kernel in the sense of GCM, ($(-r) I H I $(r))/($(-T) | $(t)) @. It is the difference between bra and ket 
that makes barrier tunneling possible and allows for the conservation of energy Eq.Q. The energy overlap kernel 
reduces to r) | H \ $(t)) owing to the choice of the overlap value ($(— r) | <I>(t)) = 1 that follows from 
{(j)k{-T) I 0i(T)) = 6ki- However, the overlap of the normalized bra and ket, (($(-t) | $(-t))($(t) | $(r)))~^/^, is 
smaller than 1. Bounce may be thought of as one of many trial tunneling paths {4'k}^ each given as two sets of wave 
bmctions, {0ifc(T)} and {(t)2k{T)}, defined on the interval [0,T/2], and related to the variables of Eqs.Q: 

, f N / 0ifc(-T) for r< 0,1 , s 

^>^^^^-{ cb,k{r) forr>0 |- 

At r = both $1 and $2 are equal to some constrained HF (CHF) state $(0) at the outer slope of the barrier with 
the constraint — 9t$(0). The Eqs.([3]) and ([5]), rewritten in terms of (pik{T) and (l>2k{T) are 

ndrCb2k + CHt) ~ ek)(f>2k = 0, (10) 

-ndr(f>ik + {h'^ (t) - ek)(l)ik = 0, 



with ft,(r) = h[(pik, 4'2k]- It should be clear that one can equally well use fields restricted to [—T/2, 0]. 

The paper is organized as follows: The main results for the HF instanton method are contained in sections IV- VII. 
These are: the variational principle (sect. IV), the formulation in terms of coordinates and momenta and comparison 
to the cranking method (sect. V), the introduction of special variables in the form of the time-even density matrix and 
the time-odd hcrmitean operator that make plain the adiabatic limit of the theory (sect. VI) and the demonstration 
that the GCM+GOA action follows from that for instanton after neglecting the velocity-momentum relations (sect. 
VII). In section VIII these results are generalized to systems with pairing. Section HI prepares useful formulas for 
later sections. Section II introduces some unusual features of the instanton method. Conclusions are given in section 
IX. 



II. GENERAL OVERVIEW 



A few comments on several unusual features of the instanton equations may be helpful. 

As the linear combination of the s.p. wave functions changes their Slater determinant only up to a factor, one 
may expect that the instanton equations may be more general than Eqs.Q which fix in a specific way lengths and 
angles among each of the sets {(pik} and {02fc} separately. Such a more general equation will imply a more general 
expression for the instanton action than Eq.([7]), and both will be given in the next section. 

The non-local in time form of the instanton equations follows directly from the transformation of the standard 
variational principle of the TDHF theory, S /(*(i) | ihdt - H \ *(t))di = Q, to imaginary time, 5 /($(-t) | hdr + H \ 
^{T))dT = 0. The Eqs.(I21) and ([5]), without the periodicity-fixing terms, have the canonical form in strange variables 

dr 64>k{T)'' 



none of which has a determined time parity. The usual canonical variables are r-even coordinates and r-odd momenta. 
Such standard coordinates and momenta may be introduced by a change of variables with the resulting equations 
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of motion local in time and canonical in form (section V). It should be stressed that a local form of the instanton 
equations does not facilitate their solution, but makes easier their comparison to other theories of the large amplitude 
collective motion (LACM). One possibility is given by [3]: (t^k — exp(— whith pk time-even and Xk time-odd. 
For one real- valued wave function and potential energy being a functional of density, V[p\, one obtains the continuity 
and "fluid velocity" equations, as for the density-phase representation of the Schrodinger equation. Energy H becomes: 



n = — dx 



V[p], (12) 



2 8p 

where the minus sign shows the role of x in the lowering of energy down to Egs in the barrier region. From the boundary 
conditions, symmetries of p and Xj and the continuity equation, one obtains action: S = (h^/m) J drdxpi^xY' ■ For 
simple systems, like the Bose-Einstein condensate of ^Li atoms, this framework allows for the exact treatment of the 
collapse of the metastable state Q- However, the density-phase variables seem unsuitable for fermions due to the 
spinor structure and the rearrangement of nodes of s.p. wave functions along the barrier that makes phases singular. 
More appropriate variables are defined in sections V and VI. 

The other peculiarity of Eqs.([3]) is that they may be thought of as describing a forced motion: The mean field h 
that causes the evolution of $(r) depends on $(— r) so one may say that one state drags the other. More specifically, 
as climbing the barrier is impossible without an external drive, the drag is necessary at the beginning of the motion 
from the metastable state through the barrier and at the beginning of the return motion to the metastable minimum. 
Action Eq. ([7]) is given by the integral of the scalar product between the change in the driven state and the state that 
drives it. For motions for which the result of the dragging is fixed by the instanton boundary conditions there must 
be some minimal dragging that causes this (fixed) result. Hence, one can expect that there is a minimum principle 
which selects instantons. If so, then solving Eqs.Q and finding decay exponent could be done by a minimization 
of a functional. The functional is practically given by Eq.Q. What remains to be done is to learn the necessary 
additional constraints which make this action minimal for instantons. 

In fission studies, mean-field states are parametrized by expectation values of observables that provide coordinates 
along the barrier, called deformations. Consider as an example the quadrupole moment Q. For bounce states <I'(t) 
one has two possible labels: Within the imaginary-time formalism, a natural choice is: Q{t) = ($(— t) | Q \ $(t)) = 
EkiM-r) I Q I Mt)). Since Q is hermitean, Q(-r) = g*(T), and Q(t) =dQ/dT = ^^((/)fc(-T) | [h{T),Q] \ M^)), 
with Q{—t) = —Q*{t). Thus the real part of Q{t) fixes r = as the return (or bounce) point. Another possibility 
is to trace deformation of the normalized state $(t), g(r) = (^(t) | Q \ $(t))/ | ^(t) p. Generally, g(r) 7^ Q{t) and 
7^ <i{~t): except for t = and ±T/2. 

Instanton cannot depend solely on a time-even variable like the real part of Q, as then Eq.Q at r = would require 
a static HF solution without constraints which cannot exist on the barrier slope. One can observe that the derivative 
q(r = 0) is equal to 2^J2kii^-r<l>k{0)]± \ Q \ </>/c(0)), where [dr4>k{0)]± is perpendicular to all 0^(0). Considering <i>(0) 
as a stationary HF state with the constraint — 9t-$(0), one can see that q{0), up to a positive constant, is the scalar 
product of two constraints: the one of instanton at r = Eqs.® and the other, — (5$(0), the proper quadrupole 
constraint on the slope where dH/dQ < 0. Since 0fc(O) lives on this slope and has the quadrupole moment (5(0), 
it must be close to some Q-constrained HF state. Hence this scalar product and the derivative q{0) are very likely 
positive. Indeed, it was found positive in the simple model 0. If so, the return point for the coordinate q{T) is at 
T > 0, which means that at t = the quadrupole moment of the normalized state $2 still increases while that of $1 
(istanton evolving backwards) decreases. Moreover, as show calculations for simple systems, states 4>ik and (/)2fc with 
the same q are different. Thus, neither Q, nor q are sufficient as labels for bounce. 

In general, the instanton mean field is not hermitean. The condition it satisfies, h[^T) = h'^{T), imposes the 
following conditions on the hermitean and antihermitean parts of its standard decomposition /i(t) — hj^^r) + hAir): 
hii{—T) = hji{T) = /i^(r), and hA{—T) = —hA{T) — h\{T). The antihermitean mean field Ha comes from r- 
odd components of densities appearing in energy ?i, either in the form of the expectation value of H or in the 
form of energy functional. In the latter case, as for the Skyrme energy functional, the generic contribution to 
Ha in the tunneling problem comes from the current density j. In the imaginary-time formalism, it takes a form: 
j(T) — 'Ek{(j)k{T)y (j)*j.{—T) — (/)^.(— T)V(/)fc(T))/2, which follows from this part of Eqs.® that shows the continuity of 
the probability flow. It follows that j(— t) = — j*(T). This differs by the factor (— i) from the conventional current in 
the real-time TDHF. As a result, the time-odd contribution to the TDHF mean field i • j • V becomes — j • V in the 
imaginary-time formalism. Its antihermitean part is proportional to the real part of j(T), and the latter appears as 
soon as the real parts of functions 4>k{T) and (/)fc(— r) become different. The time-odd mean field Ha is the immediate 
imaginary-time analogue of the Thouless-Valatin potential in TDHF [lo| . and we will use this name for it. 
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III. VARIOUS FORMS OF INSTANTON ACTION AND EQUATIONS 

The value of S which determines the fission probabihty rehes only on a part of information contained in the bounce 
solution. By using general identities: (c't-/)(— r) = — 9t-(/(— r)), f'^^dT[f{T) — /(— r)] = 0, and the constancy of 
diagonal overlaps Eq.® one can recast Eq. Q into the following forms: 

I.T/2 



S/h ^ - f drY^iMr) \ drlM-r)]) (13) 
= ^ [ drY^iM-r) \ drMr)} 

J-T/2 ^ 

= 23? / drY^iM-r) \ drMr))- 

Jo u 



The first equality shows that action for instanton evolving backwards in time, (/>fe(— r), equals to minus action for the 
instanton. The second equality shows that instanton action is a real number; the third one expresses action in terms 
of variables 0ifc and <^ik defined by Eqs.(fTO|. 

Since dr \ ^k) = {dr In | |) | + u, with v _L| (l)k{T)), and (I)i{t) for all / 7^ fc are perpendicular to <j)k{~T), the 
integrand {(f>k{—T) \ dr4>k(j)) is the sum of the full derivative plus the contribution from the component [9r0fe]_L of 
the derivative dr4>k orthogonal to the subspace spanned by all vectors {4ik{Ty}]^^i- After integration from — T/2 to 
r/2, only the latter contribution is left 

S = h drJ2{[M-r)h I [^ir)h), (14) 

J-T/2 f. 

where [4>ki^T)]± is the component of (f)k{—T) perpendicular to {4>k{T)}^^i- This shows that [0fe(— t)]_l are the essential 
variables conjugate to 4'i{t), while the components of (j)k{—T) in the subspace {4>k{T)}^^i are completely fixed by the 
overlap constraints Eq.([B]). 

As bounce ^{t) is a closed cycle in the Hilbert space ($(— r) 7^ ^(t) unlike for a line segment), action 5* may be 
written in a form of the contour integral: 

S^h<fj2{M-r)\dMr)), (15) 

which manifests reparametrization invariance of S: it does not depend at all on the instanton "speed". As can be 
seen from Eq. ljlSp . the only important features are: the path traced by | (j)k) in the vector space of s.p. states and the 
rule which associates pairs ((/)fe(— t) | and | (pkij)). Reparametrizations of imaginary time, t — > 9[t), that are both 
invertible {dr/de > 0) and consistent with the association rule: r(-6l) = -t(6I), (r(-e/2) = -T/2, t{Q/2) = T/2), 
leave S invariant. However, the reparametrized bounce, 0^(0) is not a solution to Eqs.Q. Instead, it solves: 



+ (^^ j {K9) - e,)Md) = 0. (16) 

One can recover action if bounce is known up to an r-dependent invertible linear transformation. Consider states 
^fe(T) related to bounce <pk{T) by means of such a transformation N{t): 

0fc(r)=^iV,fc(T)V'i(T). (17) 
I 

Assume N{t) = / at r = ±T/2 and r = 0. Suppose that the overlaps {ipki—T) \ ^i{t)) are given by the matrix 
M(t): 

Mkiir) = {M-r) I V'i(T)), (18) 
so that M{~t) — M{t)^ . The condition {(j)k{~T) \ (f>i{T)) — Ski means that 

N+{-t)M{t)N{t) = I, (19) 



dr 
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which leads to M ^{t) = N{t)N~^ {—t). Calculate action in terms of states ipk{T)- The integrand is: 

J2m-r){M-r)\dANu{T)Mr)]) = (20) 

ikl 



kl 



The second term is just: TrN ^drN — (3r(ln det A^). From Eq. (fT5|) one obtains: 

I.T/2 



S/h = 23? / dr V A/;-i(T)(^fc(-T) I drMr)), 
kl 



(21) 



where the omitted residual term, 3? ^'^!^^^dTdTih\detN{Ty) is identically zero, and the integration interval may be 
reduced to [0, T /2] due to the properties of M(t). Expanding either dr^'kiT') or ^pk{—T) onto ipiir) and a component 
perpendicular to all {'0fc(''")}fc=ij ons can notice that only the part [dTipk{T)]i_ orthogonal to all i^iir) contributes to 

action: S/h = 2di drj:,, M-\T){[M-r)h I [drMr)]±)- 

The Slater determinants | ^(t)), built out of ipk{T), are related to bounce determinant states | ^(t)) by: | ^(t)) = 
det A^(t) I ^'(t)), so that (*(-r) | 5'(t)) = det M(t) and H = ($(-t) | H \ $(t)) = («'(-t) | H \ *(T))/(*(-r) | 
^(t)). Therefore, energy overlap kernel TL, like action, does not involve N{t) alone and may be expressed as }11| : 

n^Y.(^,{-T)\i\i;',{T)) + ^Y.(^,{-T)^,{-r)\v\^i{T)^^^^^ (22) 

where the states iP'{t) are related to ^/j{t) via: V'K''') — 12k^ki^i^)''^k{T). The s.p. Hamiltonian h may be expressed 
in terms of various densities which do not involve N{t) either, as for example p{t) — 4'i{~'^)'4''i{T) = J2k I i^k{T) P 
+ Y.kiMki^i^M*ii-^)\^MT): etc. The Eqs. @ do involve N{t): 



i^i = 0, (23) 



but become independent of it when projected onto a space orthogonal to all {V'fe(''")}feLi' 

(ndrMr) + h{T)Mr))^ = 0, (24) 

and only this part is relevant for action. 

When the transformation N{t) has the property of a "generalized unitarity", N~^{t) = iV+(— r), the overlaps 
of states ipki—T) and V';(''") have canonical form A1~^{t) — N{t)N~^{—t) = I. Then each of the matrices NeN~^ 
and {drN)N^^ has a hermitean component which is r-even and an antihermitean component which is r-odd. For 
an arbitrary nonsingular N(t), in particular, such that keeps states ipk{T) orthonormal, the matrix M(t) in general 
depends on r and has no defined r-parity. Conversly, a general form of the instanton equation: 

ndrMr) + kr)Mr) +Y.^ik{T)Mr) = 0, (25) 

preserves overlaps Eq.® if £{t) has a hermitean r-even and antihermitean r-odd parts. There is a great variety of 
possible instanton representations with different overlaps {ipk{T) \ iPi{t)) corresponding to different matrices £. The 

periodicity condition for instanton imposes integral conditions: /5t'/2 dTdr{{'4'k{T) \ ipiiT))) = 0, i.e. integral relations 

between the matrix elements of hji, £ and the overlaps {ipk{T) \ ^->i{t)): 

( 2(V^fc(r) I hn{T) I Mr)) + Y.^£*mk{r) {Mr) \ Mr)) + (Mr) \ Mir))£mi{r))] = 0. (26) 

To assure orthonormal {ipk} at r = 0, both sets of integrals, /q^^^ and /^j^ya' should be zero. From Eqs.® wc know 
that Ski = —^k^ki provides one of the possible choices, but obviously there are many others, among them those with 
the diagonal matrix 5, i.e. with some r-dependent s.p. energies efc(r). 
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For representations with orthonormal s.p. states tpk{T), like for the usual HF determinants, Eqs.® do not hold, 
while the following relations are satisfied: ipk{-T) = J2i ^kii'^)'^ii'^) + ['^ki-T)]±, and {[ipki-T)]± \ [ipii-T)]±) = 
6ki — {M{T)M^{T))ki. Among them exists a special representation for which {ij^kiT) \ dri^i{T)) = 0, which means that 
dr as an operator has only particle-hole (p-h) matrix elements. This corresponds to the matrix £ which fulfils the 

equality £ = —h on the subspace spanned by {'4'kiT)}i^^i. 

IV. VARIATIONAL PRINCIPLE FOR BOUNCE ACTION 

Consider variation of action S in terms of some trial fission path defined in terms of s.p. states 0i/c(t) and (j)2k{T) 
for < r < T/2 as in Eq.®, fulfilling instanton boundary conditions 

5{Slh) = J2 HHki-r) I drMr)) - {dAM~r)] \ <5</>fc(r))) (27) 

k '^~T/2 

( Wifc W I 3r02fc(r)) - (a.0ife(T) I 502fc(r))) + c.c. 

k -^0 

If the states (j)2k fulfil the first set of Eqs. (fTO|l with (pit taken as the bra, then 

f.T/2 



^^^Y. [ {(^Mr) I ek - Kt) I 02fc(T)) - (hdrMr) I ^02fc(T))) + c.c. 



(28) 



k 

If, additionally, energy is kept constant so that variations fulfil: 



6H^J2 (WifcW I kr)cl>2k{T)) + {h{-T)^ik{T) I Scj,2k{r))) - 0, (29) 



then, since ((/>ifc(r) | (j)2k{T)) — 1, variation of S reads: 

SS = y^[ ' {{h{~-T) - tk)4>lk{T) - hdrMr) I S(b2k{T)) + C.C. (30) 

k Jo 

As may be seen from this equation, after the first set of Eqs. (fTO| and the condition Ti. — Egs are fulfilled, action S 
ceases to be a functional of both and (j)2k and becomes a functional of 4>2k and their time derivatives dr4'2k. The 
functions provide, through the s.p. Hamiltonian, the drive for (j)2k which is exactly required to produce dr<t>2k. 
As we have argued in section II, and as follows from the physical meaning of action, for such a driven motion S must 
be positive. Since 5S[4)2k] Eq. ((30|) vanishes for (pik that fulfil the second set of instanton equations (fTO| . i.e. when 
(j)ik and (j)2k together form instanton, the instanton action must be a minimum of 5[(/)2fc]- Thus, for (j)2k and such 
that: both fulfil the instanton boundary conditions, the overlap condition Eq.®, the energy condition Ti = Egs and 
(pik solve the first set of Egs. fTO)) for dr4>2k - calculated action provides an upper bound for action of the optimal 
(i.e. the one with the smallest action, if there are a few) instanton. That the last condition is necessary may be 
seen from the negative sign of action for bounces evolving backwards in r, Eq. (fT3|) . The assumption that half of the 
bounce equations are fulfilled eliminates trial paths with admixtures of instantons evolving backwards which would 
leave the sign of action undecided. (Note, that action for instanton evolving backwards attains the maximal among 
negative values.) This is in complete analogy to mechanics, where the real motion (g^, qi) minimizes action J Pidqi 
under the condition of constant energy provided that canonical relations Qi — dH/dpi are satisfied on each path. The 
variables introduced in the next section will make this analogy even closer. 

One can use the principle of minimal action for any representation of a trial path. A simple choice is to take for ip2k 
some orthonormal HE states with the proper boundary conditions and to look for such ["(Aifeli that ipik = 'ip2k + [V'ife]± 
fulfil Eqs. ([M)) with some r-reparametrization like in Egs. pB]) 

{ddgi:2k + Hi^ik, V'2fc]V'2fe) ^ = 0. (31) 



In this representation the overlap conditions are automatically fulfilled. Leaving r-reparametrization free one gains a 
parameter 6 that allows to control boimce velocity, i.e. the energy condition. One can decompose the s.p. mean-field 
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Hamiltonian as suggested by the formula for density p preceding Eqs. ([23)) . h[^ik,4'2k\ = 'i[^2fe] + AV'[[V'ifc]±, '02fc], 
with V the s.p. potential, so that the equation for ["01/0] ± becomes 

±,^2k]lp2k^ ^- (32) 

For complex wave functions, Eqs. ([32|) should be solved together with their complex conjugate for both [i/'ifcj^ and 
[V'ife]-L- For small [V'lfelij one could expand the r.h.s. of this equation to linear terms in particle-hole components 
[V'1-l]p/i with respect to {ip2h}, [''Pih]± — Sp[V'i±]p?i \ P >i and try to solve the system of linear equations with the 

matrix: d[AVip2h]p/d[iplj_]p'h' ■ This matrix, {pp' \ v \ hh'), where tilde means antisymmetrization, is the off-diagonal 
block of the RPA matrix (with respect to the HF state built of {ip2h}), which also appears in the ATDHF, cf Eqs. 
(2.25-2.29) and (8.24) in [§], also The solution of Eqs. (021) should be obtained for many velocities 6 to find the 
one which matches the energy condition. For larger barriers, larger differences between ipik and ip2k are necessary to 
lower the energy overlap kernel H to Egg. Then, the solution to Eqs. ((3T|) or ((32)l beyond the linear limit does not 
seem trivial. However, if found by any means, it provides action S being an upper bound for the decay exponent. 



V. INSTANTONS IN COORDINATE-MOMENTUM VARIABLES 



There are natural choices of instanton variables that correspond to time-even coordinates and time-odd momenta. 
One possibility is given by ^2i]: Mr) = ^Pkir) - ^^(t), 4>k{-T) = ((ffc(r) -I- ^^(t). It follows that ipk{-T) = (fkir) and 
^fc(-T) = -^fe(T). Due to the boundary conditions, ^fc(±T/2) = Vf^, ^^(0) = 0^(0), U±T/2) = ^(0) = 0. Thus, 
ipk are average tunneling states (coordinates) which may be parametrized by some deformation Q{t) (or its real part, 
cf section II), so that dr^k = Qdq^pk- The r-odd components ^fc must be proportional to r-odd derivative Q(t), i.e. 
to collective velocity. These two sets of states fulfil the system of equations: 

where we have used decomposition /i(t) — hn + Ha- These equations may be obtained either by decomposing Eqs. (3) 
or by deriving equations of motion from the functional J (ir(<I>(— r) | hdr + H \ ^(t)) expressed by tpk and ^k- In the 
latter case, one has to remember that ipk{—T) and £,k{—T) no longer exist as independent variables. The canonical 
form of Eqs. ([55]). without the periodicity- fixing terms, is 

n?fM . « , ,34) 

Mr) _ 5H 



dT s^HtY 

with canonical pairs {^Pk,£X) and (<^^,C/£)- Densities may be expressed in terms of ^pk and S,k^ for example, one has: 
p{x) — X^fedV'fcC^)!^ ^ ~ 2i3((p^.(a;)^fc(a;))), etc. The conserved overlaps in terms of the amplitudes ^pk and ^, 



k 



read: 



(</5fc I "^i) - {ik I 6) = ^kh (35) 
{Vk I 6) - I Vi) = 0. 



The first set of Eqs. lfSS]) is consistent with being proportional to the collective velocity Q{t). In particular, Ha 
contains S,k in odd orders, for example, the antihermitean component of the part (— j • V) of the Skyrme-type s.p. mean 
field is proportional to a piece —{{^l^Pi — iy9*V^i)/2 -I- c.c.) of the current density j. The adiabatic hmit corresponds 
to small Q and thus small \ \- 

It may be seen that the instanton dependence on Q allows to satisfy the bounce condition at r = 0: as = QCfc 
with ffe T-even, the time derivative in the second set of Eqs. ((33)) . Q^k + Q^dg^k reduces to QCfe at r = 0, where Q = 
(we assume real Q). Then hA{0) = 0, so from the first set of Eqs. ([33]) . ffe(O) — (A(0) — €k)~^dQ^pk{0)- Substituting 
this to the second set, we obtain the bounce condition at r = 



^7^(0) = (MO)-efe)Vfe(O), 



(36) 
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where Q — ^dQ^/dQ is negative at r = 0, and is determined as a function of Q by the energy condition 
T^i'PkiQ), Q^kiQ)] = Egg. The exact Eqs.([36|) foUow from the combined Eqs. ([33|) and therefore should not be imposed 
on trial paths in a variational search for instantons. 

Due to the symmetry properties of the amplitudes, action reads 



J-T/2 V \ 



d(pk \ 
dr I " 



(37) 



In this expression one immediately recognizes the familiar form J pidqi. The first set of Eqs. (|33[) are the velocity- 
momentum relations which should be fulfilled on trial trajectories in a search for bounce as a minimum of the action 
functional. Solving formally for momenta and substituting into action one obtains 



S^2h 



T/2 



T/2 



dr ( + hA{T)ipk 



k 



hnir) - ek 




(38) 



Let us compare this formula to a standard treatment of the spontaneous fission, in which one uses a family of static 
HF states, each constrained to have a prescribed quadrupole moment q, with values of q covering the barrier region. In 
such a study, one has to assume some form of the mass parameter M{q) that allows to express collective kinetic energy 
as iM(g)q^ and action as / M{q)qdq, with the implicitly understood energy conservation V{q) — Eg g. = ^M{q)q^. 

In the cranking approximation, M{q) — 2h^ J2ki'^''Pk/dq \ {had{q) — ek{q))~^ \ dipk/dq), with the adiabatic mean-field 
Hamiltonian had and its eigenenergies ek depending on q. After introducing a reparametrization q{t) in terms of some 
'time' variable t to have the correspondence with Eq. (l38p . action in the cranking approximation can be written as 



(39) 




had{t) - ek{t) 




One can see that Eq. f38)) . after neglecting the Thouless-Valatin term, is deceptively similar to the cranking expression. 
(The Thouless-Valatin term changes cranking masses by less than 20% [l^l-) However, a closer look reveals important 
diff'erences: The constants in the denominator in psp are the s.p. energies at the metastable HF minimum, not 
the adiabatic eigenenergies ek{q{t)); the states (pk, generally not orthonormal, are not equal to the adiabatic s.p. 
eigenstates ipk{q{t)); the self-consistent s.p. Hamiltonian in the instanton method depends on r-odd amplitudes, 
h — h[ipk,S,k], and this requires an iterative solution of the velocity-momentum relations. 

As follows from section HI, could be replaced in the instanton Eqs. (|10l33p by some r-dependent quantities £&(''■)• 
Such a change results from scaling the s.p. bounce states via ^^(t) = (/)'j,(t) exp(J(j (ej. — efc(r'))(iT'/?i), with r-even 
ik- This is a particular linear transformation of the type (|17p which preserves canonical overlaps Eq.([5]) and the 

periodicity, if the conditions J^^^ dTAek{T)/h — are satisfied with Ae^ — ek — ^k{T). After such transformation, 
£_k ~ cosh(j/)^^. — sinh{y)ip'f. with j/(t) = dr'Aet/h, so both and y have to be of the order Q to keep ^k ^ Q for 
small Q. This requires that the average Ae^ be of the order Q^, so only a mild deformation-dependence of adiabatic 
energies is compatible with bounce properties. 

A trial fission path is adiabatic if {(pk} differ only a little from orthonormal eigenstates of with energies efc(T) 
obtained by such a rescaling, and the velocity- momentum relations produce small (^k- Then hii[ipk] may be considered 
the adiabatic mean field and the cranking amplitudes ^k solve the second set of Eqs. (l35]) . This suggests (and will be 
shown by a different method in the next section) that in the adibatic limit Scrank provides an upper bound of (j38p 
with the neglected Thouless-Valatin term. 

Otherwise, when the larger ^k are required, the self-consistency and conditions (j35p induce a large difference between 
the contents of the cranking and instanton-motivated forms of action. For ^k not small, the enforcement of the velocity- 
momentum conditions together with Eqs. l|35p seems difficult. The same difficulty remains in the action minimization 
within this representation: since the properties of solutions to (j33p are not assured for trial paths, the conditions for 
overlaps (|35p should be imposed on them independently of other necessary conditions. 



VI. ADIABATIC LIMIT OF THE INSTANTON METHOD 



A framework analogous to that of the ATDHF theory may be obtained by defining other variables. One can observe 
that, due to the overlap conditions ([5]), a linear transformation that maps each ^^(t) into (j)k{—T) may be completed 
to a hermitean operator. Denoting the square root of this operator at each r as exp(S'(r)), with S{t) hermitean, we 
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have exp{2S{T))(f>k{T) = 4>k{—T) for all r and k. Substituting — r for r in this relation and comparing both, we infer 
that S{—t) = —S{t). Then, exp(5'(r))(/)fc(T) = exp(S {—T))4>k{—T) for all r and k. This means that the above defined 
vectors, which we will call ipokir), are time-even and orthonormal. Thus we have 

Mr) = exp(-^(r))Vofc(r), (40) 
M-t) = exp(S'(T))'(/'ofc(T), 



with "ipokiT) some r-even orthonormal states and S{t) a r-odd operator. The relation of these new variables to those 
from the previous section is given by: (pk = cosh(S')'0ofc and — smh{S)'ipok ■ The condition 5*+ = S ensures the 
constant overlaps Eq.(l6]). The bounce boundary conditions in terms of the new coordinates read: ?/'ofc(±r/2) = ipff^ , 
i^okiO) — </'fc(0) and S{±T/2) — 5(0) — 0. The states ipok define a r-even density matrix analogous to the pQ of the 
ATDHF theory [l^l- However, the object poe^ does not define any density matrix, contrary to e^^p^e'^^ of the 
ATDHF. The r-odd matrix S must be proportional to Q(r). It introduces time-odd components to the s.p. wave 
functions and its smallness is equivalent to the adiabaticity condition. The instanton equations may be written as: 

?i(e^(a,e-^)V'ofc + dr^ok) + e^(^(r) - eu)e-^^ok = 0. (41) 
Using expansions (with any operator O): 

e^Oe-^ = + [S,0] + ^[S, [S,0]] + ^[S, [S, [S,0]]] + ... (42) 

e^idre-^) = -(drS+^[S,drS] + ^[S,[S,drS]] + 



one can split Ea. (|^T|l into r-even and r-odd parts. So obtained equations are exact when the full expansion is kept. 
Since h — ft,[e"^^fcOj e~'^'(/'feo]j ^^{t) contains all even, and hA{T) all odd orders of S. The approximation valid to the 
n-th order in S consists in keeping the appropriate number of terms in both hji and Ha in each term of the equations. 

In the adiabatic limit one expects that the time derivative introduces one order of smallness, so, for example, drS 
is of the order of 5^. Then, up to the terms of the second order in S the equations read 

[hR~ek-WrS+\[S,[SM] + [S,hA]^i^ok = 0, (43) 
fidr^Qk + ( [S, ho] -f Ka] V'ofe = 0, 



with the first order Ha, and hn of the order zero, equal to = /i['(/'ofe]i except for the first term of the first equation, 
where the second order hft, should be used. In the time-odd equation, the lacking terms start at the order three, and 
would include —^[S,drS]il)Qk^ etc. As discussed in the previous section, the difference between constants and the 
adiabatic energies £&('''), which may be understood as the expectation values {ipok \ ho \ ipok)-, resides in the diagonal 
part of drS, generically of the order Q^. Clearly, not every static HF path is a proper candidate for r-even bounce 
components ipok, even if bounce is adiabatic (i.e. S is small). 
In terms of V'ofc and S action is given by 

S/h = V(V'ofe I e^idre-^) \ ^ofc), (44) 

as the part of the integrand involving dr'ipQk is identically zero due to the normalization of V'ofc- 

The approximation analogous to the ATDHF consists in solving the second Eg. pS)) up to the first order in S. With 
a given Hamiltonian, energy up to the second order in S reads Ho + 5(^0 I [S, [S, H]] \ ^0): with Ho = (^'o I H \ '^o)- 
The term quadratic in S is negative and equal to Tr(po[S, [S, Hq] + ft.A])/2, with ho = h[po] and hA linear in S. The 
latter operator is defined through its matrix elements between arbitrary states | a) and | /3): 

{a\hA\f3)=J2 (^HSM I V I /3^k) - {aiJok \ v \ p{S^Jok))^ , (45) 
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with tilde denoting antisymmetrization. Up to the second order in S, action is given by S/h = —3? X5r/2 Sfe(V'ofc I 
drS I ipok), which may be expressed as 

S/h = 23? / ^ V(Vofc I S I dr^ok). (46) 

J-T/2 

The lacking terms start at the order four, as the contribution of the order three, with the time-odd integrand 
—Tr{pQ[S,drS])/2, vanishes. After using the second Eg.ipS)). action in the adiabatic limit reads 

S = -l dTj2{{i'ok\SihA + [S,ho])\^Pok)+c.c}j ^~ f dTj2{iJok\[S,hA + [S,ho]]\iJok), (47) 

and hence is equal to the integral of: -(*o I [S, [S, H]] \ ^q) = -2{H - Ho). 

If one has an energy functional instead of Hamiltonian, one still obtains action (|Tf|) . The integrand may be shown 

equal to ~2{H[(pk,£,k] - '^[V'ofe]), with (pk = (1 + S'^/2)'i/;ofc and (,k = Sipok'- One calculates SH = Hli^ok + Sipk,£,k + 
S^k] — ^["00;!;,^*;] foi' ^Vk = 5'^'0ofc/2, ^k = S-ipok and 6£,k smaller than ^fe, to the second order in S by using Eqs. 
(El,®: 

Sn = J2 I hAir) I fk) - {S^k I hoir) \ Cfe) + {S^k \ h{r) \ fk) + cc.) . (48) 

fc 

Then one deduces S{Y,k{S.k \ fiA \ (fik) + cc.) = 2{Y,k{S£.k \ \ fk) + cc.) and S{^k \ \ ^k) = {{5ik \ K \ ^k) + cc) 
at the second order in 5*. Thus, either with the Hamiltonian or the density functional, one obtains the same form of 
the positive integrand, which, when presented as x mass, defines a positive mass for tunneling. 

In ATDHF, the linear response limit of the time-odd equation, i.e. the counterpart of the second Eq.(f4l. is: 
hdr'ipok + {ihi + [x, ^oDV'ofc — 0, with Hq — h[ipok], hi — iTr2(w[x, Po]), v the antisymmetrized interaction and Tr2 
meaning trace over coordinates of the second particle. However, hi = —ihA{x), so that the r-odd equation for the 
instanton operator 5 is a copy of the ATDHF equation, with S — x. Thus, in the adiabatic limit, instanton action 
defines the ATDHF mass hTr{Spa)/Q^. In both cases only the particle-hole components of S are determined. 

The first Eq. (|43|) provides the adiabaticity condition for a trial path, as in ATDHF [1, [l3] , but with a different 
sign by drS. It is worth emphasizing though, that this condition was practically never checked in calculations of 
ATDHF masses. Thus, up to now, decay probabilities were calculated without knowing whether a chosen fission path 
is compatible with this equation. As far as action is concerned, the difference between the real- and imaginary-time 
dynamics, i.e. between oscillations and tunneling, appears in the next order. 

A search for instanton in the adiabatic limit would consist in looking for the minimum of action determined by 
the ATDHF mass over trial paths that should fulfil the adiabaticity condition. It is well known that near the s.p. 
level crossing at the Fermi surface, an extremely small velocity is needed to keep the occupation of the lower level. 
Since in ATDHF Q must be also adjusted to keep the bulk energy Ti. constant, it may fail to fulfil two requirements 
simultaneously in the vicinity of the crossing. Thus, the proper ATDHF fission path should avoid such crossings. 
Fission paths that break many symmetries, along which crossings are avoided by a strong interaction between levels, 
could provide one remedy for this problem (as suggested by the calculations reported in 0]). The other would be 
to solve Eq. (|43p for instanton to the higher order in S, which would modify the mean field hn and avoid crossings 
present for the initial fiQ. Finally, a partial remedy is given by pairing. 

VII. BOUNCE ACTION VS GCM INERTIA 

The use of the variational principle for instantons depends on the ability to impose the velocity-momentum condi- 
tions. These conditions are crucial, as without them action for a trial path may be lower than that for bounce. Below, 
we show that the GCM formula for a collective mass that restricts generating states to r-even Slater determinants 
respects only the energy condition and hence is incompatible with the instanton method. 

Consider a family of orthonormal states labeled by the quadrupole moments gi(r) and 92 (t), t > 0, and calculate 
action Eq. (|2ip . Through the barrier, q2{T) must be different from qi{T) to make energy overlap kernel (^'(gi(T)) | 
H I 5'((72(''')))/(5'(<Zi (t)) I 5'(q2(T))) equal to Egg. If we suppose that "if depends solely on q and not on q, as in 
many GCM studies, the matrix M{t) becomes a function of qi and q2, the integrand in Eq. ((2T|) becomes equal to 



12 



Tr{M{qi,q2)-\dM{qi,q2))/dq2) and 

From this equation one can deduce a connection between the signs of S and q2 — qi- The Eqs.Q and pO| tell that 
the deformation gi of the state 4'(qi) drags deformation q2 of 5" ((72)1 thus q2 lags behind gi on the way from behind 
the barrier to the metastable minimum, i.e. 92 (t) > <1i{t)- Therefore, increasing q2 while keeping qi fixed increases 
separation between qi and 52, and thus decreases the overlap detM(qi,52)- Hence, the integrand in is negative, 
as is differential dq2 (as q{0) > q{T/2)), so action S is positive. 

In the above reasoning we used the property of the bounce equation. While using variational principle, one might 
exchange the states ^1 and ^2, and then, by the previous reasoning, negative action would follow. One might try to 
take I S* I for action in such a case, and there are cases in which this way of proceeding defines a minimum. At the 
same time, it is clear that some additional conditions are necessary in the variational formulation. 

One can expand the integrand in Ea. (|49p with respect to the quadrupole moment difference s = q2{T) — (Ji(t) 
around the midpoint q = {qi + q2) /2. When one assumes the GOA: In dot M{qi, (72) ~ ^7(9)5^/2, and then disregards 
quadratic and higher order terms in s, one obtains: 

S^-2h dq2l{q){q2-qi), (50) 

Jq(0) 

where, as discussed above, q2 > 91(92), and 7(g) = J^ki'^q'^k I dgipk) - Y,ki'\^<i'^k \ '4'i){il'i \ dqipk)- The integration 
variable dq2 = dq + ds/2 may be changed to dq, as the integral sds = d{s'^)/2 between the endpoints with s — 
vanishes. The difference of the quadrupole moments may be calculated from the constraint on the energy overlap 
kernel: Egs — H[q2,qi] ~ q] — s^{Hxy — 'Hxx)/^, where we have used the symbolic notation for derivatives of 
e.g. Tixx — d^H{x,y) \x=y=q, etc., and conditions Tix = Tiy, Hxx — Ti-yy holding for time-even H (cf. [l3|, where 
the discussion of those is given). Since the diagonal value of the energy overlap is just "potential energy" V{q) in the 
standard approach, we obtain: 



«(o) 



27(9) 



S^2h dqJ2{V{q)-Egs)[ ^ ''^'^ , (51) 

Jq(T/2) V yrLxy—rLxx/ 

where the quantity in the second parenthesis under the square root sign is the GCM+GOA mass (cf. [13] )■ 

Since additional constraints can only increase the minimum of a functional, the GGM mass must produce smaller 
action, and thus smaller decay exponent, than that of instanton. Any other action obtained with additional constraints 
will also produce larger decay exponent. As the ATDHF respects the velocity- momentum conditions to the same order 
to which it is exact, it will produce larger S than GGM. The results of calculations seem to support this, see e.g. 
[isl [l6j . On the other hand, it is known that by introducing velocities (or momenta) as additional generating 
coordinates, one can show the equivalence of such a more general GCM and the ATDHF [sl. [l7|. 



VIII. INCLUSION OF PAIRING IN THE INSTANTON METHOD 



It is well known that pairing interaction should be taken into account if realistic estimates for fission probabilities 
are to be found. In fact, it is pairing that gives the main contribution to the mass parameters, as it couples s.p. levels 
of different symmetries when they cross at the Fermi level. At the same time, it produces the gap in the quasi-particle 
spectrum which makes the collective motion more adiabatic. The proper self-consistent formalism to include pairing 
in the instanton approach is the HFB theory, in which the Slater determinants are replaced by the quasi-particle 
vacua, the many-particle states of undetermined particle number, annihilated by a set of operators: 

"^ = E(^M^°^' + ^>^)' (52) 

where operators refer to some fixed s.p. basis. We give here elements of the instanton method for systems with 
pairing. These include the imaginary-time version of the TDHFB equations, the counterpart of the formula Eq. (j2ip 
for action in terms of familiar HFB states, equations in coordinate- momentum variables (analogous to Eqs. (|33p ) and 
the formulation in terms of a time-even generalized density matrix and a time-odd hermitean opoerator that leads 
naturally to the adiabatic limit. 

For our purpose it is helpful to notice that the above customary definition implies that the HFB vacuum | 4') ^ 
exp(i Z^^a+a+) | 0), with Z = B*A*^^, depends on matrices A* and B* , while (^I' |, the corresponding bra, 
depends on A and B. 
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A. Imaginary-time TDHFB equations 

The TDHFB theory is buih on the condition of unitarity of the time-dependent Bogolyubov transformation and 
the variational principle. The HFB transformation for imaginary time, t — > ~iT, becomes: 



a'[T} \ _ ( A^{t), B^{t) \ I a 



a{-T) ; -\B+{-t), A+{-t) 



a 



(53) 



where A(t) and B{t) became functions of r, while their complex conjugate A*{t) and B*{t) became functions of 
— T. The unitarity of the HFB transformation in the real-time formalism translates to the following condition in the 
imaginary-time version: 



B+i-T), A+{-t) ; - I B*{-r), A{t) 



(54) 



This equation means that fermionic anticommutation relations for operators a+,ai, transfer to: {ai{—T),aj{—T)} — 
{af {t) , aj' (t)} = 0, and {ai(— r), a^(T)} = 6ij (and vice versa). Denoting A/'(t) the imaginary-time HFB transfor- 
mation Ea. ([53)) . its properties may be concisely written as A/'~^(r) = A/'+(— t) = a^M^ {t)(Jx, using the Pauli matrix 
notation for the block matrix. Written as separate conditions these are eight matrix equations which reduce to four 
independent relations in which r may be both positive or negative: 

A+{-t)A{t)+B+{~t)B{t) = /, (55) 

A^[t)B{t)+B'^{t)A{t) = 0, 

A*{-t)A^{t)+B{t)B+{-t) = /, 

A*{t)B'^{-t)+B{-t)A+{t) ^ 0. 

The first of those differs from the usual HFB condition as it forces anticommutation between annihilation and creation 
operators of two different sets of r and — r. This means that the usual relations: {^^(t), a~''(T)} = 5ij are not ensured. 
However, as shown below, new operators related to a(±T) may be defined, fulfilling usual conditions. 

The variational principle that gives TDHFB equations, transformed to imaginary time t — > —it, becomes: 
5 J dri^ir) \ hd/dr + H \ 'i>(-r)) = 0. Calculating variations 5/5A*^^{-T) and 6/6B'^Yt) one has to use, as in the 
real-time case, the transformation conditions Eqs. ((55|) and account for the resulting redundancy of the variables A 
and B. The term with the time derivative that defines action becomes: 

S/h^ J dT($(r) \dr^i-T)) = \J dTTr[drA+{-T)A{r) + drB+{~r)B{r)] (56) 

= -\j dTTr[A+{-T)drA{T)+B+{-T)drB{T)]. 

The matrix element of Hamiltonian ($(t) | H \ <&(— r)) is expressed by contractions: 

mr)\a+a^m-T)) = p^..{r) = {B*{-t)B^ {r))^,, (57) 
($(r) I a.a,, \ <i>(-r)) = k^,{t) = {B* {~t)A^ (r))^,, 
($(r) I a>+ I <i>(-r)) ^ ^^.(r) ^ (A*(-r)B^(r))^., 

which, due to conditions (j55p . have the following properties when regarded as matrices: 

p(-r) = p+(r), (58) 

Using those and proceeding as in the case of TDHFB we arrive at imaginary-time TDHFB equations written symbol- 
ically (where only the second index of the amplitudes is explicit): 

Ak(r)\ i + t(r). A(t) \/ Au(r)\ A,(r) 
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where, for a given Hamiltonian, the self-consistent potential: r^^(r) = '^^givp^-y^s — v^j,~fSv)ps-y{T) and the pairing 

potential: /S.^u{t) — 'Yl-^s'^t^'^iS'^i^i'^) have the properties: T{—t) = r+(r), and A"^(r) — — A(t). The same properties 
hold for the mean fields with additional rearrangement terms that follow from a density functional. These ensure the 
property h{-~T) ~ /i^(t) of the mean-field Hamiltonian /i(t) = t + r(T), and the same property, h(— r) — h^(T) of 
the total HFB mean-field Hamiltonian h(T) given by the matrix in Eqs. ((59| . As a result of this, the equations ((59|) 
conserve both energy and all relations (|55p . The terms with constants Ek on the r.h.s. fix the periodicity of solutions 
and these constants are equal to the quasi-particle energies at the metastable HFB ground-state. The bounce solution 
to Eqs. ([59|) has to be periodic and provide a path connecting the HFB ground state | gs) with some HFB state 
I $(t = 0)) at the same energy beyond the barrier. 



B. Variational principle 



In a similar way as in the HF case, one can deduce the minimum principle for action under conditions of constant 
energy and fulfilled Eas. (|59|) for < r < T/2. The redundancy of variables A, B complicates the Hamilton equations, 
but the following relations hold: -25n = Efc(('5Wfc(-r) | Mr) \ Wfc(r)) -I- (Wfc(-T) | h(r) | 5Wk{T))) and -255 = 
I drWkir)) ~ (a^[Wfc(-r)] | 5Wfc(T))), with Wfe denoting the vector composed of Since 
taking a formal variation of S* + 7Y with respect to (5>V^ and SV^k leads to the correct equations (|59l) , the arguments 
of sect. IV can be repeated and one obtains the same constraints that specify bounce as the minimum of action (note 
(>Vfc(-T) I W,(t)) 

The first of Eqs. ((55)) means that ($(r) | $(— t)) = 1. Since these two HFB states are different, the imaginary-time 
HFB transformation determined by the matrices A(±t) and B{±t) cannot be unitary. However, it may be related 
to a normal unitary HFB transformation given by some matrices U{t), V{t) via some invertible, though non-unitary 
matrices C(t). Let us suppose a relation: 



c^tir) = }_^C,.{T)f3+{T), (60) 



3 

with quasi-particle creation operators Pf {t) related via some U{t) and V{t) matrices to a+, a^t, namely [cf Eq. (|53p 



a+(r)\ / (t/(r)C(r)f, (nr)C(r)f \ ( a+ 
c^i-r) j [ (T/(-r)C(-r))+, (t/(-r)C(-r))+ i I a 



(61) 



It follows that U{t),V{t) define the same Z{t) as A{t) and B{t) do and that U{t)+U{t) + V{t)+V{t) = 
C+^^{t){A+{t)A{t) + B+{t)B{t))C~^{t). If one chooses C(t) that transforms the hermitean matrix A+{t)A{t) + 
B^{t)B{t) to the unit matrix, then U{t) and V{t) become matrices of a standard HFB transformation. Now, the 
first of Eas. (|55)) means that: 

{U{~t)+U{t) + V{~t)+V{t))-^ = C(t)C{-t)+, (62) 

while three other follow from this and from the HFB properties of matrices U (r), V{t) and U{—t), V{—t). The second 
equation ([55| is just the condition of the antisymmetry of Z{t), the equations three and four: (/ -|- Z'^{t)Z{—t))~^ + 
Z+{t){I + Z{-t)Z+{t))''^Z{-t) = I and the antisymmetry of matrices: Z{t)+{I + Z{-t)Z+ {t))'^ and (/ + 
Z{—t)Z^{t))^^Z{—t), follow from the previous two. 

Using the same reasoning as the one leading to Eq. (PT|) . instanton action ([SS]) can be expressed in terms of the 
normalized HFB states | ^(r)), defined by U{t) and V{t), using relation ((62)l : 



1, 



2 



T/2 



S/n = --^ / dTTr[{U+{^r)UiT)+V+i-T)ViT))-\U+i-T)drUiT) + V+i-T)drViT))], (63) 



T/2 



where we have omitted the integral of dr Indet C(t) between the endpoints, as it is purely imaginary. 

The contractions Eq. ((58|l can be expressed through C/(±t), V{±t) and the corresponding HFB states \E'(±t) in the 
following way: 

p,. = iV*i-r)mrV)-^V-ir)U = l\r" , (64) 



1 at 


1 




-r)) 






-r)) 




1 


flp 1 


*(- 


r)) 


(*(r) 




-r)) 
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-r)) 



(*(r) I ^i-r)) 



15 



where the matrix U{t) = {—t)U (t) + V~^{—t)V{t) is related to the overlap of standard HFB states via: (^'(t) | 
^(-t)) = [detC/(r)]i/2 [li]. 

Now, one can treat ([63|) as a functional on trial fission paths ^'(t), defined by two families of HFB states ^i(t) 
and *2(r) for < t < T/2 

^ ' ^ [ *2(t) forr >0, / ^^^^ 

smoothly connecting some HFB state $(0) beyond the barrier at energy Egg to the metastable ground state '^gs, and 
fulfilling the condition of constant energy overlap and the Eas. ([59|) for5'2(r). Taking [/(r) = {t)U2{t)+V^ {t)V2{t) 
for T > 0, and having U{t) = [/+(— t) for r < 0, one can calculate action as: 

S/h = ~^ [ ' dTTr{l]-^{r){JJt{T)drU2{T) + V^{T)drV2[T))\. (66) 
Jo 

The minimization of this action over fission paths that fulfil constraints should reproduce the bounce action. Its value 
for a trial path that satisfies constraints is an upper bound for the bounce decay exponent. 



C. Coordinate and momentum variables 



The coordinate-momentum variables may be introduced in a similar way as in section V. Decomposing amplitudes 



into T-even and r-odd components, A(t) = A+(r) — A- 



B{- 



B+{t) + i?-(T), with and B+ matching 



(r), A{-t) = A+{t) + A_(r), B{t) = B+{r) - B^{t), 
5 at T = ±r/2 and $(0) at r = 0, and A_ = B_ = at 



T = 0, ±r/2, one obtains the system of equations (with only the second index of the amplitudes made explicit) 




hR - Ek, 



-A_, 
A+, 



hR — Ek, 
-hA, 



A 



-h 



-h\-Eu, -A! 



+ 

R-Ek 











) 









(67) 



with the mean fields h — Hr + Ha and A = A+ + A_, with A+(— r) — A+(t) and A_(— r) — — A_(t). In a similar 



way as for Eqs. (|33p. the first two Eqs. (l67l) connect velocities drA^k, drB+k with momenta A^k and B^k, showing 
that they all, together with the r-odd mean-field potentials Ha and A_, are proportional to the collective velocity Q. 
In the coordinate-momentum representation, these are the constraints that must be imposed on trial fission paths to 
assure that bounce provides the minimum of the action functional. The Eqs. ([55]) written in terms of new amplitudes 
become eight relations which may be combined to four r-even and four r-odd equations, e.g. the first Eq. l(55|) leads 
to A+A+ - AtA_ + B+B+ - B+B^ = I and AtA+ - A+A_ + B+B+ - B+B^ = 0, etc. 

Let us call the diagonal and off-diagonal submatrices of the matrix in Eq. (j67p — and h.R. From symmetries and 
definitions it is clear that hij(r) is hermitean and time-even and h^(r) - antihcrmitean and time-odd. In imaginary- 
time TDHFB, the operator is the generalization of the Thouless-Valatin mean field Ha of the ATDHF. 

Denote the vector built of A+k and B+k as 9^ and the one built of A^k and B^k as S^, i.e. yVfc(r) — 0fc(r) — Sfc(r). 
Then the Eqs. ((671) take the form 



UdrQk = -h^efe + (hfl - £'fc)Sfc, 
hdr^k = {iiR- Ek)Qk-i^A'E.k■ 



m 



The variation of energy written in terms of Ofc and :^k reads 

2Sn = J2 {{S^k I h^Sfc) - {SOk I hROk) - {SEk \ h^e,-) + {SEk \ h^Sfc) + c.c.) . (69) 

k 

The three last terms, together with their complex conjugate, contribute at the second order in r-odd components, 
assuming and 6E,k being of the first, and 6Qk of the second order of smallness. Owing to the r-parity of the 
amplitudes, and after integrating by parts, action reads: 

f ' dTTr[At{T)drA+{T)^Bt{T)drB^{T)\. (70) 
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This can be expressed a.s S — —U^ \ drQk), i-e. it is imaginary-time TDHFB action in the form 

J Pidqi. Substituting from the first Eq. (|68p one can obtain the cranking-hke expression for action as in sect. V. 

D. Adiabtic expansion and limit 

The above formulas are a copy of those in sections V and VI, up to the common factor (—1/2) appearing in the 
expressions for S and STi.. Hence, after showing that the operator that maps amphtudes at r onto those at — r is 
hermitean one could represent HFB bounce in terms of r-even amphtudes and a r-odd hermitean operator S, as in 
sect. VI, and repeat the whole reasoning on the adiabatic limit of the instanton method. (To emphasize the analogy, 
we keep the same notation for the time-odd operator as in HF, although it acts in the enlarged space.) 

The argument goes as follows: The HFB transformation from operators (a^(T), a(— r)) to {a'^ {~T),a{T)) is 
N{~t)M^^{t) — N{~'t)M^ {—t) (cf Eq. ([51)) ). indeed hermitean. Calling this transformation exp(25(T)), with 
S{t) hermitean, and considering its inverse, we have S{—t) = —S{t). Then, we find that exp(iS(T))A/'(T) = 
exp(5(— r))A/'(— r), so calling this r-even transformation M{t), we have J\f^^{T) = A/'+(t), so M{t) is a regular 
HFB transformation. Denoting its amplitudes u and v, we have 

-as?.)- 

The properties of N{t) and N{t) imply a^S'^ {T)ax = —S{t). As we need a relation between amplitudes and these 
form columns of the matrices A/'^(t) and Sfjij), we notice that A/"^(t) = [A/'^(t) exp(-5^(T))(A(^)"H'^)]-^^('^)j 
and that the matrix A/'^(t) exp(— 5"^)(t))(A/'^)~^(t) is hermitean owing to the HFB property of A/'(r). Moreover, 
due to this property, one has A/'^(t) = exp(-5'(T))A/'^(r) with the hermitean, r-odd S{t) = M"^ {t)S'^ {t)M* {t) . It 
follows from the properties oi N and S that axS'^ {T)ax = —S{t). Thus 

^= (4* 4*)' ^^^^ 

with s+ = s, and = — f. With this S'(t), we have the expected relations 

( b:!:;! ) - ( :{;! ) • ( ^ ) - =p<-*(^» ( ) ■ (^^) 

where only the second index of the amplitudes is shown. With these, all the results of the section VI can be repeated 
for imaginary-time TDHFB. In particular, the integrand of the action integral S*, which in terms of the amplitudes 
Wofe = (wfc,Ufc) and the operator S reads — | ^^^((^tWo*: | S \ Wofc) + c.c.) is equal to —2{H — Ho) at the second 
order in S, hence positive. The Eqs. take exactly the form (|¥T|l of the imaginary-time TDHF, with obvious 

replacements of Wok for ipok and h for h. They reduce to the form at the second order in S. 

The TDHFB equations may be also formulated in terms of the generalized density matrix. The counterpart of the 
HFB density matrix in the imaginary-time formalism is 



p(r), «(r) \_fB*{-T) 



(B^(r),A^(r)) = (a, exp(5*(T))a,)7^o(T)(a, exp(-^^(r))a,), (74) 



with 72-0 (r) the HFB density matrix corresponding to A/'(t). Owing to the property of S, it is equal to 
e-'5(^)72o(r)e'5(^). This non-hermitean quantity, call it TZ, apart from not being any HFB density matrix, is an ana- 
logue (note that TV = TZ) of the density matrix in the ATDHFB theory [ll|, 7^ = e'^TZoe-'^. Iri^ terms of it, Eqs.dSH) 
read: hdrTZ + [h,7?.] = 0. The r-odd part of this equation, linear in S, obtained by expanding TZ = TZq — [S,TZa] + ... 
and discarding the second order quantity [ho,7?.o], 

hdrTZo + [[S, ho] + iiA, TZo] = 0, (75) 

is an alternative form of the second Ea. (|^5)) in terms of 7?.o and S. Its solution is identical to the ATDHFB solution, 
S = x- This follows directly from the structure of the building blocks of the Thouless-Valatin mean field Iia. One has 
hA = Tr{vpi) and A_a^ = Y.^ys'^aH-iSKi-iS, with pi = -[s,po] " ^^o^^*, = pof + f{pl - 1) - skq - kqs* . Since, 

in ATDHFB, TZi = i[x,TZo], one has = ihi, where hi is the ATDHFB time-odd mean field for x = S. Thus, the 
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adiabatic TDHFB instanton method produces mass given by: mass x Q"^ = ^Tr{'R.oS), equal to the ATDHFB mass, 
of [l^. In the zero pairing hmit this mass reduces to the ATDHF value hTr{pos) / . 

A reasoning similar to the one presented in section VII shows that, within the GCM approach, a use of some r-evcn 
pairing variable (for example, the pairing gap) as a generator coordinate, without fulfilling the velocity-momentum 
relations, will lead to a smaller decay exponent than that for bounce. 

IX. CONCLUSIONS 

We have presented the instanton method for nuclear fission in various representations. This has allowed for some 
comparisons with other methods commonly used in fission studies. We have also sketched the imaginary-time version 
of the TDHFB theory, which allows to include pairing. 

There are many similarities between the instantons describing quantum tunneling and the periodic TDHF solu- 
tions. Both appear as a result of the quasi-classical approximation, find a natural formulation in terms of time-even 
coordinates and time-odd momenta and reduce to the same time-odd ATDHF equation in the lowest order in mo- 
menta. The ATDHF equation for a path, which should be fulfilled for consistency, is usually not checked for static 
paths constructed by means of the CHF. When the velocity- momentum equations require small momenta that violate 
energy conservation, this means that the chosen path is far from instanton. 

The main difference between the two methods is that in quantum tunneling there is no single HF state or density 
matrix, but one deals with two different states, bra and ket. This happens to be the very reason for the existence of the 
minimum principle: it defines the minimal driving of one state by the other, necessary for tunneling. Instanton action 
turns out to be a minimum of the action functional when the constraints of constant energy and velocity-momentum 
relations are imposed on trial fission paths. Action calculated for any such path would provide an upper bound for 
the decay exponent. We argue that the ATDHF (ATDHFB) mass respects those constraints, while the GCM-I-GOA 
mass does not. The main practical problem is how to construct trial paths fulfilling the constraints. 

The need for two Slater determinants for instanton leads to another important difference between the mean-field 
studies of oscillations and quantum tunneling: The instanton method relies on the off-diagonal matrix elements of 
the Hamiltonian, which are beyond the usual scope of the mean-field theory. To use instantons in practice, one has 
to define various off-diagonal matrix elements of the commonly used effective interactions, like, for example, of the 
density-dependent term of the Skyrme-like force (for its possible definitions see [1^ ) and of the Coulomb-exchange 
interaction. 

When comparing the instanton method to theories of large amplitude collective motion (LACM) one has to recognize 
that the aims of the latter are much wider than those of the former [UIMl- In LACM, equations for the collective 
path or action are a source of formulas for potential and inertia tensor of an effective Hamiltonian in a restricted set of 
deformation coordinates and conjugate momenta. Often the next step consists in the requantization. The supposed 
universality of the so conceived effective theory for LACM underlies the whole procedure. On the contrary, instanton 
should be found once for a studied decay. No interpretation of the integrand in the action formula as mass x is 
necessary. It could be even dangerous, as in some representations of instanton these integrands are piecewise negative. 
Only the value of the integral has the physical significance and this does not depend on the representation. 

Of course, one could extract collective inertia from action represented with a positive integrand, but the positivity is 
obvious only in the adiabatic limit. In a general case, action Ea. (|38p contains momenta to all even orders, and the 
higher order terms become naturally more important for higher barriers. Hence one expects that mass also depends 
on the barrier height, or energy, when tunneling from excited states is considered. A small energy dependence of mass 
is seen even for highly collective Bose- Einstein condensate Q- 

For pairing gaps of ~ 1 MeV and for not too high fission barriers, the adiabatic approximation may be satisfactory 
for many fission paths. Then it may appear that the most important in the search for instanton is the exploration 
of a sufficiently rich family of paths, preferably with as few preserved symmetries as possible, while ATDHFB action 
(including Thouless-Valatin terms) is a sufficient estimate of the instanton action. 

Even if this is true, fission of odd-Z or odd-A'^ nuclei will require much more effort to understand, within the 
instanton method, a dramatic significance of the odd fermion and of the specific mean fields induced by it that break 
time-reversal invariance. 

It is clear that the method considered here is applicable to quantum tunneling in any fermion system, provided it 
has a meaningful mean-field description. Extensions to include thermal effects and decay form excited states seem 
also straightforward. The real progress of the method will depend on practical solutions. 
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